Stata：coefplot 如何选择交互项系数估计值进行绘图？

为了让大家更好的理解本文的内容，欢迎各位培训班会员参加明晚 8 点的直播课 「Stata：coefplot 如何选择交互项系数估计值进行绘图」

最近有小伙伴问了这样的一个问题：

也就是使用 coefplot 展示系数估计值与置信区间时如何选择交互项变量。

方法一：查看所有可用的变量名称

这里我们使用 auto 数据集为例：

在回归中使用因子交乘算子：

直接使用 coefplot 就可以绘制所有变量的：

那么该如何仅仅选择交互项的系数绘制呢？

回归的结果也会被存储在返回值中，查看 e-class 返回值：

回归系数都会被存储在 e(b) 中：

提取这个矩阵的列名：

这些列名都是可以在 coefplot 的 keep() 选项中使用的：

需要注意的是，e(b) 中值为 0 的不能用，这些是被 omit 的变量：

也可以使用 foreach 循环提取 e(b) 中值不为 0 、含有 # 的进行绘图：

方法二：使用通配符

还可以使用通配符：

另外提取回归系数的时候也可以遵循这个方法流程：

所以虚拟变量的系数可以这样提取：

方法三：从返回值中提取数据使用 twoway 直接绘图

对该方法感兴趣的小伙伴也可以学习下这个课程：

「Stata：如何展示回归系数和置信区间」（https://rstata.duanshu.com/#/brief/course/2d3950fe01ed46119619b615dc8a90e4）

coefplot 虽然方便，但是有时候可能不能完全满足需要，所以也可以直接使用 twoway 绘制：

由于 eret list 中的 e(b) 中只存储了回归系数，不足以直接绘图，所以下面我们还是使用 ret list 中的 r(table) 数据绘图：

r(table) 的列名实际上也是和 e(b) 一样的。

然后我们再把这个数据转换为 dta 数据：

然后就可以使用 twoway 绘图了：

如果想给 y 轴添加标签的话，可以这样：

矩阵转换成 dta 也可以直接使用 svmat：

然后后面的操作就一样了。

方法四：使用 serset 提取图表数据修改后再重新绘制

对于这部分内容感兴趣的小伙伴还可以前往系列课程「Stata 编程导论」的第 15 课时：Do file 编程实战（六）学习，里面提供了更多的案例。

系列课程「Stata 编程导论」：https://rstata.duanshu.com/#/course/500094e7213744bc904e4e2cb270d135 Do file 编程实战（六）：https://rstata.duanshu.com/#/course/27a3b7a55a2545a0a6281ae8fafaf15d

使用 coefplot 绘图之后我们可以使用 graph save 把图表保存为 gph 文件：

重新调用：

为了下面使用方便，编写一个 setnames 命令：

gph 格式的图文件会保存绘制其的代码和数据，数据可以使用 serset 调用出来：

查看绘图代码：

稍微改动好看点：

可以看到交乘项是 __00000F = 2 或者 __00000F = 3

然后就可以重新绘制了：

这样我们就非常彻底的解决了这个问题！

直播信息

为了让大家更好的理解上面的内容，欢迎各位培训班会员参加明晚 8 点的直播课 「Stata：coefplot 如何选择交互项系数估计值进行绘图」

1. 直播地址：腾讯会议(需要报名 RStata 培训班参加)
2. 讲义材料：需要报名 RStata 培训班，详情可阅读：一起来学习 R 语言和 Stata 啦！学习过程中遇到的问题也可以随时提问！

更多关于 RStata 会员的更多信息可添加微信号 r_stata 咨询：

附件下载（点击文末的阅读原文即可跳转）：

https://rstata.duanshu.com/#/brief/course/e877d6b6020447e89530f36639d10c6f